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Abstract: A linear function submission-based double-auction (EES-DA) mechanism for a 
regional electricity network is proposed in this paper. Each agent in the network is equipped 
with a battery and a generator. Each agent simultaneously becomes a producer and consumer 
of electricity, i.e., a prosumer and trades electricity in the regional market at a variable price. 
In the EES-DA, each agent uses linear demand and supply functions when they submit bids 
and asks to an auctioneer in the regional market. The EES-DA can achieve an exact balance 
between electricity demand and supply for each time slot throughout the learning phase 
and was shown capable of solving the primal problem of maximizing the social welfare of 
the network without any central price setter, e.g., a utility or a large electricity company, 
in contrast with conventional real-time pricing (RTF). This paper presents a clarification of 
the relationship between the RTF algorithm derived on the basis of a dual decomposition 
framework and EES-DA. Specifically, we proved that the changes in the price profile of the 
EES-DA mechanism are equal to those achieved by the RTF mechanism derived from the 
dual decomposition framework except for a constant factor. 

Keywords: Multi-agent system; Distributed algorithm; Dual decomposition; Double 
auction; Eagrangian relaxation 
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NOMENCLATURE 


/;+ G oo) 

Electric energy consumption profile 

G [ 0 , 4 -’”“) 

Electric energy generation profile 

b[^ G 

Battery charge profile 

b\~ G 

Battery discharge profile 

m;+ G [0,m;+’”“] 

Profile of electric energy sold to the local electricity market 

m. G [0,m. ] 

Profile of electric energy bought from the local electricity market 

G [0,oo) 

Profile of electric energy sold to the outside grid 

gr e [0,gr-n 

Profile of electric energy bought from the outside grid 

A 

Profile of state vector 

s\ G 

Profile of the state of charge (SOC) of the battery 

li e [0,1] 

Storage efficiency 

7G [0,1] 

Electricity transmission efficiency 

c\ 

Cost function for generating electric energy 

D\ 

Utility function for consuming electric energy 

<^1 

Individual utility function 

wl 

Individual welfare function 

Pt 

Price profile 

P?^ 

Price of electricity sold to the outside grid 

P?- 

Price of electricity bought from the outside grid 

a; 

Constant term of parameters of the bidding function 


Primary coefficient term of parameters of the bidding function 


The superseript t & denotes the t-th time slot in a day. The subseript i G denotes the i-th agent 
in an eleetrieity network. 

1. INTRODUCTION In reeent years, solar power and other renewable forms of energy, ineluding 

wind and hydroeleetrie power, have been attraeting attention beeause fossil fuel will beeome diffieult to 
mine in the near future. In eontrast with thermal power generation based on fossil fuel, eleetrie power 
generated by photovoltaic (PV) cells and wind turbines is uncontrollable and unpredictable because it 
depends on the weather and climate. Therefore, adapting consumer demand to the fluctuating supply 
is a key issue when developing our future electricity network. Recently, a decentralized autonomous 
electricity network has been studied as a way to overcome this problem [1]. 

In this paper, we propose a double auction-based demand-side management mechanism for a regional 
prosumers’ electricity network named the inter-intelligent renewable energy network (i-Rene) [2]. A 
schematic view of i-Rene is shown in Figure 1 . In i-Rene, locally generated electric power is consumed 
in the local electricity network. Currently, surplus energy produced by PVs connected to the electricity 
network in Japan is sold at a fixed high price, i.e., the price is set as the Feed-in Tariff, and transferred 
to a grid managed by a regionally monopolizing electric power company, almost forcibly. This causes 
reverse power flow, which could destabilize the global electric power network. 
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Figure 1. Overview of inter-intelligent renewable energy network (i-Rene) 

Therefore, i-Rene is designed not to transfer surplus energy to the external grid, to whieh it is 
conneeted via a gateway, and which essentially transmits electric power in a single direction. Each house 
is equipped with a smart meter running intelligent software. The software automatically trades electricity 
in the regional electricity market and allocates the generated electricity appropriately by referring to the 
load profile of each house. The regional electricity market is automatically managed by the gateway. 
We assume that the regional electricity network covers an entire town in which the number of houses is 
assumed to be between ten and several thousand. 

We assume that each house connected to i-Rene contains a generator using renewable energy 
resources and a battery. The intelligent software in the smart meter automatically trades electricity 
on behalf of the dweller. However, local power generation, trading, and electricity stored in the batteries 
cannot always satisfy the total demand in the local grid. In such cases, consumers connected to i-Rene 
would have the option of buying electricity from an outside power grid through the gateway at a high 
price. i-Rene is not centrally managed in that there is no central utility or large electricity company; 
instead, the agents in the network simply trade their electricity through double auction. Therefore, 
designing a double-auction mechanism and analyzing its properties are important. 

Double-auction mechanisms have been used previously for designing multi-agent software systems 
to manage regional electricity networks [2-6]. Double auction involves the submission of bids by 
potential buyers and the simultaneous submission of asks by potential sellers to an electricity market. 
If an adequate price balancing demand and supply is determined, the auctioneer clears the market. The 
determined price allows the participants to buy or sell electricity on the basis of their submitted bids 
or asks. One example of a double auction-based multi-agent software system for regional electricity 
network is the PowerMatcher [3,4,7]. Several field tests proved the ability of the PowerMatcher to 
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alloeate energy resources efficiently in a decentralized manner. In the PowerMatcher, each agent submits 
continuous demand and supply functions directly to an auctioneer in a market. Several multi-agent 
systems for realizing intelligent electricity networks employ such a function submission-based double 
auction (FS-DA) [2-4, 7, 8]. The advantage of FS-DA is that the auctioneer determines an appropriate 
clearing price by calculating the intersection between aggregate demand and supply functions without 
iterative communication between the auctioneer and the agents. This calculation is easily performed 
using simple agent software without any manual intervention by humans. 

On the other hand, a number of studies have been devoted to the use of real-time pricing (RTP) 
in smart grids [6,9-17]. These studies mainly used game theory and/or control theory to develop 
RTP algorithms and to prove their optimality and/or robustness. In particular, the dual decomposition 
framework provides a sophisticated explanation and theoretical foundation for RTP algorithms [18]. 
However, the disadvantage of most of the dual decomposition-based RTP algorithms is that an iterative 
communication between a central utility, which settles the electricity price, and the agents are required to 
balance demand and supply. In addition, most of the dual decomposition-based RTP algorithms assume 
that a central utility “does not” behave selfishly, but behaves so as to optimize social welfare, i.e., it 
attempts to solve the master problem dedicatedly. 

The use of the double-auction mechanism leads to the automatic generation of a time-varying price 
profile by simply balancing demand and supply. The consumers and the power producers are expected 
to change their load, generation, and storage profiles by reacting to the time-varying price profile so as 
to maximize their profit. This means the mechanism has an intrinsic RTP mechanism. However, the 
optimality and the performance of LFS-DA in a regional prosumers’ electricity network with distributed 
power resources and batteries have not been proved. 

In this paper, we describe a linear function submission-based double auction (LFS-DA) mechanism 
for a day-ahead power market, in which an adequate price profile is gradually formed to control the 
demand and supply without any central price controllers. Throughout the learning phase, the LFS-DA is 
able to achieve an exact balance between electricity demand and supply for each time slot. We show the 
LFS-DA to be essentially equal to dual decomposition-based RTP from the viewpoint of changes in the 
price profile. Especially, our main contributions are as follows. 

• We propose an LFS-DA mechanism for a regional prosumers’ electricity network. The mechanism 
is able to achieve an exact balance between electricity demand and supply at each moment during 
iteration. 

• We prove that the changes in the price profile obtained by the LFS-DA become the same as those 
obtained by the RTP derived from the dual decomposition framework except for a constant factor. 

Our main result, i.e.. Theorem 1, is very simple but powerful. It shows that simple LFS-DA is 
sufficient for managing the day-ahead market in i-Rene without the need for a central utility to perform 
RTP to balance demand and supply, and attempts to maximize social welfare by setting an appropriate 
price profile. 

The remainder of this paper is organized as follows. Section 2 provides the background to our 
research and related work. Section 3 describes the problem definitions and basic assumptions of 
the target electricity network used in our research, which is a regional electricity network consisting 
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Figure 2. (Left) real-time prieing in a conventional power grid, and (right) double-auction 
pricing mechanism in a regional prosumers’ electricity network 


of many prosumers. It also describes an RTF algorithm derived from the dual decomposition 
framework. Section 4 describes our proposed LFS-DA mechanism and provides proof that the LFS-DA 
is intrinsically equivalent to dual decomposition-based RTF. Section 5 contains details of a simulation 
experiment and its results. Finally, Section 6 concludes this paper. 

2. BACKGROUND AND RELATED WORK 


2.1. Penetration of distributed energy resources and smart grids Distributed energy resources (DER), 

e.g., distributed energy generators and storage facilities, have gradually been introduced into households, 
plants, and offices. Such distributed generation (DG) is mostly performed by using renewable energy 
resources, including solar, wind, and hydroelectric power. Conventionally, especially in Japan, large 
power generation companies have generated electricity centrally with thermal electric power stations 
and nuclear power plants, which they then transmit to consumers on a unilateral basis. The situation 
differs from country to country. Specifically, in Japan, each power supplier has traditionally been a 
regional monopolizer, obliged by law to completely fulfill consumers’ electricity demands. 

However, the increasing penetration of DERs and the emergence of smart grid technologies are 
changing the requirements of power networks. The wide-scale introduction of DERs causes a reverse 
power flow to conventional power grids, which would be rendered increasingly unstable and confused. 
Therefore, it is important to develop a regional electricity network capable of mitigating the effects of a 
reverse power flow. In this context, several power network concepts, including a smart grid [19], micro 
grid [20], and a virtual power plant (VFF) [21], have been studied with the aim of smoothly introducing 
DERs to the grid. 
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In the development of a regional network based on DERs, demand side management (DSM) becomes 
an important design element. However, in an electricity network, balancing the demand and supply at 
each moment is crucially important. Therefore, achieving a balance between the demand and supply 
requires the power demand, rather than the power supply, to be controlled. A comprehensive survey of 
research relating to DSM was conducted by Alizadeh et al. [9] . Moreover, many field tests demonstrating 
the effectiveness of DSM have already been performed, e.g., [7,16,22]. 

2.2. Demand side management The realization of efficient, robust, and applicable DSM has led to the 

study of market-based methods [9], which have used price-based incentives, e.g., time-of-use (TOU) 
pricing, RTF, and peak-time rebates, to control consumers’ electricity demand. The approaches used in 
game and control theory have provided the theoretical foundations for DSM methodologies. 

A comprehensive survey of the application of game theory to smart grid systems was provided by Saad 
et al. [6]. Vytelingum et al. analyzed the Nash equilibrium for an electricity grid to which distributed 
micro-storage devices are connected. Their multi-agent simulation showed that this had the effect of 
flattening the electricity price profile [23]. Mohsenian-Rad et al. proposed a distributed algorithm for 
autonomous DSM [17]. They provide analytical proof to show that their proposed algorithm guarantees 
a significant reduction in the total energy cost of the system. 

Control theory forms the theoretical foundation of RTF on the basis of a dual decomposition 
framework [10-15]. A comprehensive tutorial of dual decomposition was compiled by Falomar et 
al. [18]. In a dual decomposition-based approach, the problem of maximizing the social welfare of 
consumers and/or a power producer under several constraints is defined as the primal problem. This 
problem is transformed into a dual problem by applying a Lagrange relaxation and decomposition 
methods which have historically been used in optimization problems in electric power systems [24]. 
This dual problem can be divided into many sub-problems and a master problem. Each sub-problem 
corresponds to a selfish profit maximization problem of a consumer or producer, whereas the master 
problem corresponds to the problem involving a search for the optimal price profile. In a dual 
decomposition approach, the introduced Lagrange multipliers are interpreted as “prices”. In this way, 
the dual decomposition method provides a mathematical foundation based on RTF and derives each 
pricing algorithm of DSM depending on each condition of the target electric power network. On 
the basis of the dual decomposition framework and related concepts, various RTF methods have been 
proposed [10-15,25-28]. 

2.3. Multi-agent system and double auction Most of the previous studies relating to RTF control theory 

assume that power producers and consumers exist separately in an electric power network as shown in 
Figure 2 (left). In contrast, in our work, we focus on an electricity network consisting of electricity 
prosumers as shown in Figure 2 (right). The term prosumer refers to a person or entity who is both a 
producer and a consumer at the same time [29]. If a consumer installs a FV system at his/her house, 
he/she not only becomes a consumer, but also a power producer. The economical structure of an energy 
network in which regional prosumers participate would differ completely from that of a conventional 
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unilateral power grid. Such a network has been studied in the context of multi-agent system-based 
microgrid. 

In recent times, multi-agent systems using the double-auction method for managing electric power 
networks have been gaining attention [20,30]. Kok et al. developed the PowerMatcher Smart Grid 
Technology [7], a multi-agent-based distributed software system for the market integration of small 
and medium sized DER units. In a PowerMatcher cluster, many agents interact with each other and 
the equilibrium price for electricity is determined automatically by using the double auction approach. 
The effectiveness of dynamic pricing using the PowerMatcher has been tested in simulations and field 
experiments [3,4,7]. Vytelingum et al. proposed a market-based mechanism on the basis of Continuous 
Double Auction (CDA), and the use of an efficient trading agent strategy. They showed that the 
mechanism and the agents were able to cope with unforeseen demand or increased supply capacity 
in real time [5]. Taniguchi et al. proposed an inter-intelligent renewable energy network (i-Rene), which 
is automatically managed by using double auction [2]. Related papers reported the capability of their 
methods to achieve a reduction in peak consumption and/or the efficient utilization of energy [2,5,7]. 
However, they were unable to prove that the double auction-based approaches maximize the social 
welfare of the prosumers’ network from an analytical point of view. This contrasts with the previously 
mentioned dual decomposition-based methods, all of which have an analytical foundation that guarantees 
an increase in the social welfare. 

Most of the multi-agent-based approaches employ double auction in the regional electric power 
network to determine the price and amount of electricity transacted. In double auction, bids (i.e., 
buy orders) and asks (i.e., sell orders) have to be matched to enable a certain amount of electricity 
to be transacted and the price to be determined. However, if various agents submit bids and asks 
at various prices for various units of electricity, the market auctioneer has to use a large number of 
iterative communications to match demand and supply.To reduce the communication overhead, Kok et 
al. simplified the communication by using only demand functions and prices [7]. Agents in the network 
submit individual supply and/or demand functions to the market and if all of the agents submit continuous 
functions, the market is able to calculate an adequate clearing price and suitable transactions by searching 
for the point at which the aggregate demand and supply functions intersect. In this paper, we refer to this 
double-auction method as the function submission-based double-auction (FS-DA) mechanism. In the 
Agent-based Modeling of Electricity Systems (AMES) test bed model developed by Ei et al. and in the 
i-Rene developed by Taniguchi et al., the submitted functions are also restricted to linear functions [2,8]. 
We refer to an FS-DA mechanism that only uses linear functions as an EFS-DA mechanism. 

In contrast with the traditional top-down approach, in which price profiles are determined by a public 
utility or by RTF generation (see Figure 2 left), our double-auction approach determines price profiles 
dynamically through interaction between sellers and buyers in a bottom-up manner (see Figure 2 right). 
Several related mechanisms were proposed. Sam et al. proposed a function bidding mechanism for 
multi-agent charging of electric vehicles [31]. Their mechanism could eliminate an iterative exchange 
of messages by introducing function bidding. Kiani et al. proposed a dynamic market mechanism 
and analyzed the market as a dynamical system [32,33]. Xu et al. proposed and analyzed a simple 
market mechanism using supply function bidding [34]. Their supply function differed from that of 
our model. The research most closely related to ours was conducted by Papavasiliou et al., who 
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Regional smart grid 


Figure 3. Power flow in the inter-intelligent renewable energy network (i-Rene) 


proposed a Newton algorithm-based double-auetion mechanism [35]. Their mechanism, supply function 
bidding, is interpreted as a Newton algorithm for optimizing problems with a decomposable structure. 
The mechanism enables an automated demand response control system through a double auction-like 
mechanism. However, they did not consider a prosumer network in which each consumer uses batteries 
and renewable energy resources, and did not demonstrate the mathematical relationship between dual 
decomposition-based RTP and the function submission-based double-auction mechanism under the 
condition that each agent behaves selfishly, i.e., they required each agent to use a Newton algorithm 
to update its variables. In this paper, we show that the price change in LFS-DA becomes the same as 
dual decomposition-based RTP except for a constant factor, which is our main theorem. 


3. PROBLEM DEFINITION 


3.1. Basic assumptions In this work, we define a target regional electricity network as follows (a 

schematic diagram is shown in Figure 1 and Figure 3): we do not distinguish between electricity suppliers 
and consumers; rather, we assume that all people on the grid are capable of producing and consuming 
electricity as prosumers. Thus, it is assumed that each house is equipped with a generator powered 
by renewable energy, e.g., PV cells or wind power generators, and a storage device, e.g., a battery. 
Furthermore, each house is considered to have a smart meter running software, which automatically 
controls the electricity of the house by managing its battery, transmitting electricity to other houses, and 
communicating with other information systems. 

In the regional electricity network, houses are connected through a regional electricity market through 
which each prosumer can sell and buy electricity. The local electricity price fluctuates on the basis of 
regional demand and supply. The prosumer attempts to optimize its trading rule, consumption, and 
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generation so as to maximize his/her welfare. Basieally, i-Rene is intended to behave as an independent 
electrieity network. 

The outside grid provides the prosumer with an optional alternative from which to buy electricity at 
the fixed high price at any time. In addition, the prosumer can also sell surplus electricity to the outside 
grid at the fixed low price, which is set so as to inhibit reverse power flow. We assume the outside grid is 
a unilateral conventional grid, as is common in Japan. Hereafter, we refer to a household, including the 
prosumers living in the house and the software, simply as an agent. 

Suppose there are N agents in i-Rene and a set of these agents is represented by ^ := {1,2,..., V}. 
An agent can consume, generate, charge, discharge, buy, and sell electricity through its smart meter 
during every one of the T time slots. In addition, the number of time slots for transactions is shared by 
all of the agents in the grid. The set of time slots is defined as ^ := {1,2,..., T}. The /-th agent can 
determine the amount it consumes /■’*", generates l\~, charges , discharges b\~ , sells to , and buys 
from the regional market , and sells to g^f^, and buys from the outside grid during each time slot 
t. The variables of the z-th agent (z G JV) at time t ^ are defined as shown under Nomenclature. The 
relationship among the variables is also schematically shown in Figure 3. 

Each variable has its individual lower limit and/or upper limit owing to the limited capability of each 
apparatus or system. The superscripts ■+ and ■ represent the direction from the viewpoint of a smart 
meter, which is at the center of electric energy flows in an agent’s house. Specifically, ■+ and ■“ represent 
outflow and inflow, respectively. We assume that the energy flow through the smart meter adheres to the 
law of the conservation of energy for each time slot t as follows. 

- r + - b\- + - m\- +g\+-g\-^Q (1) 

The amount of electricity demanded has to be balanced against the supply in the regional electricity 
network for each time slot t. We assume 

^ ( 7 zn|+ - m\-) =0 ( 2 ) 

where yG [0,1] represents the electricity transmission efficiency. If 7 = 1, there is no electricity energy 
loss during transmission. In addition, the storage efficiency rji G [0,1] must be taken into consideration. 
If rii = 1, the charged electricity can be fully charged without any loss. The storage profile s] represents 
the state of charge (SOC) of the z-th agent’s storage device at time t and is expressed by the following 
equation. 

= C+ E W 

where is the initial SOC of the z-th agent’s battery. 

Agents can buy electricity from and sell electricity to the outside grid at the price of and 
per unit, respectively. The situation for which reverse power flow is completely prohibited is modeled 
by setting pf^ = 0. An immediate resale behavior between an agent and the outside grid is suppressed 
by assuming the following constraint. 


0<P?^<P?~. 


(5) 
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3.2. Social welfare We define the cost for generating electric energy for the i-th agent by C■ : M —)■ M, 

where C- is a convex function of class and define the utility for consuming electric energy /■“ for the 
i-th agent by : M —)■ M, where D\ is a concave function of class C^. The welfare W/: x —)■ M 

of the z-th agent is defined as follows. 


Wfxup) := Y, W-iA^Pt)^ ( 6 ) 

■= 0/(4) -Ptym'i^ +Ptm\^, (7) 

0/(4) := DK/D -C/(/r) +P?+g\+-p?-g]- (8) 

■= {A)te:7, (9) 

4 := dO) 


where p = (pr)ze5^ ■— {PI:---:Pt), Pt is a buyer’s price profile in the regional electricity market at time 
t G and (^\{xt) is a utility function for the z-th agent at time t including payments to an outside grid. 
Because of the electricity losses in transmission, the seller’s price and buyer’s price become different. To 
balance the amount of money paid by sellers and buyers, the seller’s price becomes ypt. The z-th agents 
behavior at time t is represented by concisely. We call a state vector for the z-th agent at time t. 
Next, we define the social welfare of the network by W (x,p). 

ie.yVteS' 

+ Ptym\^ - Ptnfr) 

= L ( L ^K^i)+pt L (r^r-w^D) 

'-V-' 

=0 from equation (2) 

where x = (x,),g^ (xi, ...,xzv). This shows that the social welfare inside of the regional prosumers’ 
electricity network does not depend on the price profile of the internal regional market. Therefore, 
maximizing W{x,p) is the same as maximizing 0(x) := L/g^0i(x0 where 0(x) is the sum of the 
individual utility functions 0z(x,) := Ejgjr 0/(4)- 

3.3. Real-Time Pricing Algorithm The following problem is defined as the primal problem, which has 
the purpose of maximizing the social welfare of the regional network. 

Problem 1 (Primal problem). 

maximize V 0z(x/), (13) 

subject to Xi E Vz G ^, (14) 

E//W) = 0 (15) 

4 ( 4 ) ■=ym\^-m\~. 


( 11 ) 

( 12 ) 


(16) 
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Algorithm 1 Iterative solution of dual problem: real-time prieing (RTP) algorithm 
k^O 

Initialize the priee profile 

repeat 

// Eaeh agent solves its sub-problem (45), and obtains its solution. 

Update ^x*{p^^^). 

H A eentral utility solves the master problem (47) under the eondition that xi = . 

Update ^ — 0[^t{p^^^), for each t. 

k i — k \ 

until a predefined stopping criterion is satisfied, 
return Transact {x\^^}i(z^ with p^^'> as price profile 


where is a feasible set for the z-th agent except for the network constraints (see APPENDIX A). 

The primal problem is numerically solved by a solver. However, no electric power companies, 
utilities, or governments can directly control the behavior of any of the agents. In addition, they cannot 
obtain information about C- and D^- because this information is usually private. Even if the central control 
system could obtain the private information, solving the primal problem requires a huge computational 
cost when the number of agents N becomes large. In reality, solving the primal problem centrally 
is impossibly difficult in the current economic environment. Therefore, a decentralized optimization 
method using price information is gathering attention in the context of smart grid management. 

Based on the dual decomposition framework, an RTP mechanism maximizing social welfare, in which 
each agent maximizes its welfare selfishly, can be derived, and the price profile is updated centrally using 
a sub-gradient. The RTP algorithm is outlined in Algorithm 1. 

(k) 

The price profile p) is updated by using a known sub-gradient method [ 10,1 8] Adopting this method, 
Pi’^^ is updated as follows: 



(17) 


(18) 



= E (Tmr+(p'‘') -<“(?'“)), 

(19) 


where 0^ > 0 is the learning rate of the sub-gradient method. Iterative updating of Xi and p provides the 
numerical solution of the dual problem (see A). This solution is also a solution of the primal problem (13) 
- (16) . Therefore, the social welfare of the network is expected to be maximized. Eurther details are 
provided in APPENDIX A. Erom the viewpoint of the pricing mechanism, solving the dual problem 
equates to the determination of a price profile by a utility, so as to maximize the social welfare. Therefore, 
we refer to this algorithm as an RTP algorithm. 

Although this algorithm presents a feasible solution for the decentralized energy dispatch problem, 
it is unable to achieve an exact balance between demand and supply before converging to an optimal 
solution, despite the necessity of such a balance in an electricity network. In each step, 5m^ := 
(7znf+(p(^)) — mf ^ is not exactly a zero vector in most cases. Therefore, a heuristic 
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proeess should be applied when aiming to eompensate the differenee dm'. One of the simplest methods 
is that the utility eompensates this differenee when it is not zero by selling or buying dm' to or from the 
outside grid at the price of or respectively. In this work, we assume that the i-Rene gateway 
employs this heuristic procedure if i-Rene adopts the RTP algorithm. 

In practical situations, the iteration cannot be performed many times, and only a small number of 
iterations are permitted. In such cases, the effect of dm' cannot be ignored. In this context, a mechanism 
which would be able to achieve an exact balance between the demand and supply is desirable. The 
LFS-DA is such a mechanism as it has the ability to balance demand and supply exactly for each 
time slot and every iteration step, and is guaranteed to solve the problem in the same way as the dual 
decomposition-based RTP algorithm. 

4. THE LINEAR FUNCTION SUBMISSION-BASED DOUBLE AUCTION The LFS-DA 

mechanism works as follows. A prosumer submits their individual linear demand and supply functions 
parametrized by two parameters, i.e., aj and /3/, described in (20) and (21). Then, the market-clearing 
price pt is determined exactly by calculating the point at which the aggregate demand and supply 
functions intersect. After several iterations between an auctioneer and the prosumers, all prosumers 
transact their electricity using the determined clearing price. In this section, we formulate the LFS-DA 
mechanism and disclose its theoretical properties. 

4.J. Transaction with LFS-DA In LFS-DA, each agent has a linear demand function {pt) and 

a linear supply function m'^^ = ii'j^{pt) for each time slot t. These two functions are determined by using 
the two parameters a\ and /3/. 

m;+=At;+(p,):=LA'A-«/J, (20) 

= K~iPt) ■= [-PiPt + Otj \, ( 21 ) 

where [x\ := max(x,0). In the LFS-DA, each agent submits the two parameters, i.e., a' and /3/, to the 
auctioneer in the regional market to provide information about the linear demand and supply functions. 
In this work, we assume /3/ is a fixed positive constant, and a' is a flexible variable that is optimized by 
each agent. 

When each agent submits the two parameters (aj, /3/) to the market, an adequate price pt is determined 
that would suffice to clear the market by searching for the point at which the aggregate demand and 
supply functions intersect. The constraint for balancing demand and supply (2) becomes 

E fiixhPt) = 7 L Lfe - a/J - L L-A'a + 

= 0 . ( 22 ) 


CC- 

The price pt divides the agents into two groups. If for the i-th agent, -^ < pt, the agent becomes 

^ od 

a seller, i.e., Pjpt — ccj = m\^ > 0 and m\~ = 0. In contrast, if for the z-th agent > pt, the agent 
becomes a buyer, i.e., m'^^ = 0 and —ji'pt + ocj = m'-~ > 0. We define the set of suppliers at time t by 
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{Pt) '■= I ^ < PtJ £ ^}, and the set of eonsumers at time t by l^{pt) '■= I ^ > PtJ G ^}. In 
this ease, on the basis of the eonstraint (22), an appropriate elearing priee has to satisfy the following 
equation. 


( 22 ) ^0 = 7 ^ (Plpt-aj)- ^ {-pjpt + a^) 

By solving the above equation for pt, we obtain 


(23) ^pt = 




+ a[ 


'It {.Pt) 


(23) 


(24) 


E(p.) +(?.)“■'■ “f (K) ■= Em ■= Etp.) 

The priee pt satisfying (24) exaetly balanees demand and supply, and fulfills the eonstraint in (22). 


Lemma 1 (The elearing priee). The clearing price pt that satisfies (24) uniquely exists and be calculated 
exactly. 

Proof. flPiPi ~ OJ/J ~ l^PiPt + Otj\ is a pieeewise linear and monotonically inereasing function of pt 
and its value at p, = 0 is negative. Therefore Wi Pt — E/g^ L“ Pt + J ^ piecewise 

linear and monotonically increasing function of pt and its value at = 0 is negative. Therefore, the 
solution of (22) uniquely exists and can be calculated exactly. □ 


In the LFS-DA, when the auctioneer in the market receives (aJ,/3/),-e^ ,gj^ from all of the agents 
at the k-th iteration, the auctioneer determines the price by solving (24). Based on the price 

each agent either transmits = p-’*'(pf*''’^^) or receives = p-“(pf^'*’^^). After 

determining the amount of ), each agent has to reconfigure its for the given 

) so as to satisfy the constraint x,- e Under the condition that each agent behaves 
rationally, i.e., selfishly, the process of the LFS-DA mechanism is described as Algorithm 2. 


4.2. LFS-DA as an bottom-up RTP method In the same way as RTF, communications between an 

auctioneer and prosumers are required by the LFS-DA to maximize social welfare. When a price profile 
{Pt)teA" is observed, each rational agent tries to maximize its welfare by optimizing aj on the basis of 
the announced price profile (pfteA'^- 

When /3/ and pt are fixed, a} can be determined uniquely from (m'f^,m\~) as follows. 

= PiPt + ■ (25) 

Equations (20), (21), and (25) show that aj and {m\^pn\~) have bijective relation under the condition 
that /3/ and pt are fixed. Therefore, the optimal solution af(p\^^) can be obtained using the following 
procedure. 

A* ) = argmax ( (xi) + ^ pf ^ (7mJ+ - m \^) | (26) 

(pf^(pf^)). (27) 


1 


Note that /3/ is fixed in this work. 
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Algorithm 2 Iterative update in LFS-DA 
k^O 

Initialize the priee profile 

repeat 

// Eaeh agent solves its sub-problem (26), and obtains its solution. 

Updates* = 

Update af ^ -f {mf- - mf+) . 

H Eaeh agent submits (a*^/3/) to the market. 

Update ^ market_clearing((al,/3/),-e^,fe^). 

Update {pi^+^))) 

H Reeonfiguration by eaeh agent 

Update xf ^^^ ^ proj {x* \ m\^ = mf^,m\~ = mf ^^), for eaeh i 
k i — k \ 

until a predefined stopping criterion is satisfied, 
return Transact with p^^^ as price profile 


where xf = (/f +, lf ~, bf +, 




*t+ *t~ 

56 ; 56 ; 


andx* = . Equation (26) corresponds 


to an optimal solution of the sub-problem of the dual problem introduced in the APPENDIX A. 

Here, we formalize each agent’s iterative learning process. The change in the price profile and agents’ 
behavior are determined as shown in Algorithm 2. In the k-th iteration, the price profile is assumed to 
bep(*) := Each agent selfishly maximizes its welfare and obtains x*, before obtaining af = 

The clearing price settled by market_clearing can be obtained by solving equation (24). 

On the basis of the new price profile the amount of electricity for transacting through the 

regional electricity market {m^-^,m\~) is updated as follows. 


/ (/:+!)/+ (/:+!)/- 

m; ,m; 


4- 




(28) 


After the update, the state vector becomes 






,m; 


,m. 


(fc+l)r- *t+ *t- 

T Si •) Si 


(29) 


This state vector usually violates the constraints, i.e., xf ^ Therefore, each agent has to reconfigure 
to ensure satisfaction of the constraints under the condition that ). 

The map proj^.(-|m-^ = ) denotes a projection to ^i satisfying {m\^,m\~) = 

) for all t ^ 2^. In the paper, the reconfiguration is performed by solving the 
sub-problem with the additional constraints, i.e., (m|^,m|~) = ). 

The following Eemma 2 indicates that the changes in the price profile obtained by market_clearing 
can be regarded an update formula as follows. 
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Lemma 2 (Change of price profile in LFS-DA). Price profile of the LFS-DA is updated as follows. 



(30) 

&(/>“’) := E PWriPl) 

(31) 

ieJC 



(32) 


where 0^ = (/JJ 

[Pt 




(Pt. 


Proof. For all t & 


(fe+i) 
Pt 


yoc 




' Eipr^’) 


“h CC 










a, 


*t 


+'^i^ir(pr"A 








(33) 

(34) 



(35) 


+ 

+ 

1 

+ 

77-iH' irL 

E <-) 

(36) 

ieJ/' 

(37) 



□ 

Theorem 1. This update formula (30) of the LFS-DA is the same as that of the dual decomposition-based 
RTP algorithm (17) except for the constant factor 0|. 

Proof. Comparing (17) with (30) confirms the statement to be true. □ 

Therefore, if we consider the determination of the price in LFS-DA as a type of price update (17) 
in dual decomposition-based RTP, the price change in LFS-DA becomes equal to that of RTP where 

In summary, in LFS-DA, if the auctioneer simply settles a clearing price profile to balance demand 
and supply, it implicitly solves the master problem of the dual problem (47) in the same way as in the 
RTP algorithm. This means that the changes in the price profile obtained by the LFS-DA mechanism 
become the same as those obtained by the RTP derived from the dual decomposition framework. Note 
that the sub-problems in RTP and in LFS-DA are completely the same and are solved in the same way. 
(see (45) and (26)). Therefore, both of the master problems and the sub-problems are the same and are 
solved in the same way in the RTP and the LFS-DA except for the constant factor of the sub-gradient in 
the master problem. Our main theorem therefore indicates that LFS-DA can increase social welfare in 
the same way as the conventional RTP without requiring any central price controller. 
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5. EXPERIMENT This section describes our numerical example of the LFS-DA mechanism. In 

addition, we compare its results with those obtained for the dual decomposition-based RTP algorithm 
through a simulation experiment. 


5.1. Experimental conditions In this experiment, we set the number of prosumers to V = 20 and the 


number of time slots to T = 24. The utility functions are set to be 


a,(I) := 


nt 

cojl - 

Ml! 

261 


0<l<^ 

- - 9j 

(^<1 

\e 


(38) 


where Co\ and 6- > 0 are given constants [10,36]. In addition, we set 6- = 30 and coj = 10 for all 
t e £i",i e 

We considered all generators to be PV systems; therefore, no variable costs were taken into account. 
The generation costs C- are defined by := 0 (V/-“ G [0,/(~’™“]). We used the PV generation 

profiles that were measured at 20 houses in Higashi Ohmi city during Autumn 2010 as the maximum 
value of the per-slot energy production profile [ 37 ]. Fig. 4 shows for each prosumer i. 

We set = 0 , 5 ™“ = = 1,^?^’’™“ = 1, and r/,- = 0.7 for all i G We set the efficiency 

of transmitting electricity within i-Rene to 7 = 0.8, and the parameters of the sub-gradient method to 
6k = 6.1. We set the maximum values to = 5, and m7’™“ = 5 for all i G t & 57. We set the 

prices of the outside grid to = 20 and pf^ = 0 for all t G 57 . 

In this experiment, two baseline methods were compared with the LFS-DA. The first, RTP (dual 
composition), employed Algorithm 1 . However, the demand and supply were usually not balanced before 
the algorithm converged to the optimal solution, in which case the utility conducting the RTP would 
compensate the difference by buying the deficit electricity from or selling the surplus electricity to the 
outside grid, respectively. The cost for the compensation was assumed to be redistributed equally among 
the agents. RTP without compensation shows the social welfare before the difference is compensated. 
Although dual decomposition-based RTP could increase prosumers’ welfare, it was reduced because of 
compensation by the utility. The second method is without trading, that is, electricity trading does not 
occur in the regional market. Instead, agents attempt to increase their welfare by using their batteries and 
elasticity of demand. The optimal is numerically obtained by solving the primal problem using a solver. 


5.2. Result 

The value of the social welfare is shown in Figure 5 as a function of time. This clearly shows that 
the LFS-DA outperforms the RTP algorithm based on dual decomposition. Especially, in the early 
stages of the iteration, the LFS-DA was found to outperform the RTP algorithm significantly. Although 
the LFS-DA mechanism was able to obtain an exact balance between demand and supply, the RTP 
algorithm was unable to modify prosumers’ bids and asks reactively. The cost associated with the lack 
of compensation had a negative impact on the overall social welfare, although each agent would still be 
able to increase its own welfare. 
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Figure 4. Amount of power generated l\ by 20 prosumers’ PV eells 


The ehanges in priee profiles are shown in Figure 6. These figures show that the RTP algorithm 
and the LFS-DA obtained the same price profile except for numerical errors as our theory suggested. 
Figure 7 shows the rate at which each agent’s welfare increased compared with the without trading 
condition. As shown, LFS-DA increased the welfare of all the agents in contrast with the RTP algorithm, 
which caused several agents to reduce their welfare because of the compensation for the imbalance in 
demand and supply. Figure 8 shows the electricity consumption of three representative agents. For the 
without trading condition, the electricity price for agents is fixed at pf~. In contrast, RTP and LFS-DA 
resulted in fluctuating price profiles and leveled consumption profiles. This figure shows that LFS-DA 
functions in the same way as the RTP by controlling electricity consumption using price signals. 


6. CONCLUSION This paper presents a linear function submission-based double-auction (LFS-DA) 

algorithm for a regional prosumers’ electricity network named i-Rene. This mechanism can achieve an 
exact balance between electricity demand and supply for each time slot throughout the learning phase. It 
was proved that LFS-DA can solve the primal problem of maximizing the social welfare of the regional 
electricity network in the same way as the RTP algorithm, but without any central price setter, i.e., a 
utility or a large electricity company. This means that the price controller for RTP formed in a bottom-up 
manner as a result of the use of the double-auction mechanism can maximize the social welfare of the 
regional prosumers’ electricity network. This paper presents a clarification of the relationship between 
the RTP algorithm derived on the basis of a dual decomposition framework. 

The relationship between the dual decomposition-based RTP algorithm and the simple LFS-DA 
algorithm for i-Rene was clarified. Specifically, we proved that the changes in the price profile with 
the LFS-DA mechanism are equal to those achieved by RTP controlled by the algorithm derived from 
the dual decomposition framework except for a constant factor. A simulation experiment was used to 
demonstrate the proposed mechanism numerically. 
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Figure 6. Priee profiles determined by iteration with LFS-DA and with RTP 

Our main result, i.e.. Theorem 1, is very simple but powerful. We used it to show that the simple 
LFS-DA is sufficient for achieving day-ahead optimization and planning for i-Rene. We also showed 
that there is no need for a central utility to perform RTP to balance demand and supply; instead, social 
welfare was maximized by setting an appropriate price profile. Thus, the mechanism functions in a fully 
decentralized manner. 

In this paper, we fixed {Pj}ie^,te3^ which are the coefficients of the first order of the linear demand 
and supply functions. The parameters originally represented the price elasticity of prosumers’ demand 
and supply. Interestingly, our main result showed these parameters to perform the function of the learning 
rate in RTP. This means that the price elasticity of prosumers’ demand and supply affects the convergence 
of the iterative calculation of the price profile in i-Rene. 

Our main theorem showed the theoretical connection between research pertaining to a 
multi-agent-based microgrid using linear function-bidding, e.g., PowerMatcher [7], and a control 
theoretical approach to RTP. The connection is expected to lead to renewed collaboration and fruitful 
integration between the two research areas. This formulation also provides us with the possibility to 
extend the LFS-DA mechanism beyond its current capabilities. For example, Yo et al. introduced 
chance constraints to dual decomposition-based RTP [26] and a similar extension should be possible 
for the LFS-DA approach. In future, we aim to develop a double-auction mechanism, which could 
automatically control a regional electricity network of prosumers even in a noisy environment. 
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Figure 7. Rate of increase of each agent’s welfare compared with the without trading 
condition 


APPENDIX 


A. Dual Decomposition and Real-Time Pricing In this section, we derive the dual 


decomposition-based RTP algorithm on the basis of Problem 1, i.e., the primal problem. The complete 
description of the i-th agent’s feasible set except for the network constraint (15), becomes the 
following. 




(39) 

h'^{xi) < 0 

Vt e ^,J E {1 ,..., 16}, 

(40) 

h^^^(xi)=0 

Vt E 

(41) 


h^^(xi] 

. _ ^/+,mm 

-/'+ 

‘'1 1 

h'^{xi] 

■■= -r, 

h'^{xi) 

:= 

b'^^ixi) 



h^^{xi] 

:= -m^+, 

h'\xi) 

:= — 

h''^{xi) 

:= 


h^\xi] 

— -8\~. 






f—.max 




'J+ i,-l-,max 


h‘^\xi) :=b\--b 


',max 


i i 


h''-(xi) := 


t+ +,max 

m ■ — m- , 


h^^^{xi) := 


t— —.max 

m- — m- , 




h^^^(xi) : 


= -Si 


„mit 


k=l 

h^^^ixi) :=ll+-ll-+b^+-b\~ + m] 


t+ 


■m\ +g;+ 


■s\ , 


h^^\xi) 

fiixi) 


k=i 


k—\ _ max 

‘J 7 1 


:= ym- 


t-\- 


■m\ . 


In the above formulas, are constraints for the domains of the state vector x\. The two 

constraints for the battery profiles, and represent the storage capacity constraints. The constraint 
represents the balance of the electricity flow measured by the i-th agent’s smart meter. 
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Figure 8. Electricity consumption profile by three houses 

Network optimization problems with a decomposable structure, such as the primal problem, can be 
solved by using the dual decomposition technique [18]. This technique decomposes a primal problem 
into independent sub-problems representing agents’ selfish optimization behavior, and a master problem 
that optimizes prices to balance the demand and supply. Each of the decomposed sub-problems can be 
solved independently. The solution of the primal problem and that of the dual problem are known to 
be the same if Slater’s theorem is satisfied [38]. In addition to that, the Eagrange multipliers, which 
appear through dual decomposition, can be interpreted as a “price” in our economic system. Within the 
context of the studies of smart grids and power systems, the dual decomposition technique provides the 
mathematical basis of RTP. The dual problem corresponding to the primal problem becomes 

Problem 2 (Dual problem). 


minimize g(A), 

(42) 


s{^)-= E supL,A(x,-), 

(43) 


■= <l>i{xi) + E 

(44) 




Generally, it is known that the Eagrange multiplier A := represents the “price” of goods, 

which is constrained by an equation relaxed by A/ [18]. A comparison between (44) and (7) provides 
us with a clear understanding of this relationship. In this problem, Af represents a buyer’s price of 
electricity traded in the local electricity market at time t.The simultaneous solution of the sub-problems 
and the master problem results in the solution of the dual problem and this is achieved in practice by 
adopting an iterative optimization technique. 
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If the Lagrange multipliers A are given, the objeetive funetion (43) of the dual problem ean be divided 
into N objeetive funetions. Each objective function for each agent can be solved independently. 

Problem 3 (Sub-problems). 


maximize (x,), (45) 

subject to Xi e 

This problem is referred to as a sub-problem of the dual problem (42). From the viewpoint of game 
theory or micro-economics, if each agent can be regarded as rational, i.e., selfish, each sub-problem is 
expected to be solved autonomously by each agent. If the price of electricity in the market is given as A, 
one of the optimal solutions x*(A) to the problem (45) can be obtained as follows. 


= argmaxL,-;t(^<') 

XiE^i 


(46) 


where = {l*^+,bf+ 




The dual problem is solved by optimizing the price profile A G MJ' . For RTP formulation purposes, 
we assume there is a central utility determining the price profile in the regional electricity market. On the 
basis of each agent’s optimal strategy (x* (A))/^^ , the central utility has to solve the following problem. 


Problem 4 (Master problem). 


minimize g (A), (47) 

g(A)= ^La(4(A)). (48) 

This problem is referred to as the master problem and its solution corresponds to searching for an 
adequate price that would balance demand and supply. 

One of the most common procedures for solving the dual problem is as shown in Algorithm 1 . The 
Fagrange multiplier A/^^ is updated by using a known sub-gradient method [10,18] If we adopt this 
sub-gradient method. A/ is updated as follows in the dual problem: 


= A,® - SiKA,!*') 

(49) 


(50) 


(51) 


where > 0 is the learning rate of the sub-gradient method. By updating Xi and A iteratively, the 
numerical solution of the dual problem can be obtained. The solution of the dual problem is also a 
solution of the primal problem. Therefore, the social welfare of the network is expected to be maximized. 
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